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ABSTRACT 

We use the recent measurement of the velocity dispersion of star-forming, outer-disk knots by 
Herbert-Fort et al. in the nearly face-on galaxy NGC 628, in combination with other data from the 
literature, to execute a straightforward test of gravity at low accelerations. Specifically, the rotation 
curve at large radius sets the degree of non-standard acceleration and then the predicted scaleheight of 
the knots at that radius provides the test of the scenario. For our demonstration, we presume that the 
Ha knots, which are young (age < 10 Myr), are distributed like the gas from which they have recently 
formed and find a marginal (> 97% confidence) discrepancy with a modified gravity scenario given the 
current data. More interestingly, we demonstrate that there is no inherent limitation that prevents 
such a test from reaching possible discrimination at the > 4u level with a reasonable investment of 
observational resources. 
Subject headings: gravitation — galaxies: kinematics and dynamics 



1. INTRODUCTION 



The application of the Newtonian law of gravity 
to the dynamics of galaxies and galaxy clusters leads 
to the well-known discrepancy between the inferred 
dynamical m a ss and the d i rectly observed ma s s (cf . 
Zwickvl Il93l ISmithl IJ9p IRoestad fe Shostald IJ971 



Rubin. Ford, fc Thonnardlll978t ). This discrepancy is, in 
turn, the motivation for the dark matter hypothesis — 
that the vast majority of matter in the Universe is of 
some undetermined composition. Although the need for 
some dark matter in certain environments is no w well es- 
tablished (|Clowe et al.ll2006t I Angus et al.ll2007l ). and ex- 
pe riments appear tantal izingly close to directly detecting 
it (jBernabei et alJl2~008ft . the possibility that gravity de- 
viates from Newtonian behavior in the low acceleration 
regime is not precluded by the existence of dark matter. 
Much of the work on proposed alternatives to New- 
tonian gravity focus on demonstrating the consistency 
of such models with observations, particularly with the 
measured kinemati cs of disk galaxies at large galac- 
tocentric radii (cf. iBeeeman. Broeils. fc Sanders! 119911 : 
iMcGaugh fc de Blok 1998). The ultimate aim however 
is to discriminate between the modified gravity and dark 
matter hypotheses, and this requires identifying a regime 
in which the two hypotheses make distinctly different 
predictions. There are some basic requirements for what- 
ever environment one envisions constructing such a test 
in. First, one will need to work in a regime well beyond 
the scale where Newtonian gravity first appears to fail. 
Second, one will need to work in a regime appropriate for 
the proposed modified gravity model. Because the cur- 
rently proposed modifications to Newtonian gravity are 
not yet complete theories, there is little to be gained by 
invalidating a specific gravitational variant beyond its in- 
tended working parameters. Third, one will need to iden- 
tify an environment in which the physical distributions 
of the baryons and hypothesized dark matter are suffi- 
ciently different that they create distinct gravitational 
signatures. Lastly, one will need to identify a dynami- 
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cally simple environment where the interpretation does 
not depend on one's understanding of a complex dynam- 
ical history. 

All of these requirements drive us to study the out- 
skirts of disk galaxies. The observed outer-disk ro- 
tation curves provide one of the principal arguments 
for deviations from Newtonian gravity (if one assumes 
no dark matter, e.g., modified Newtonian dynamics or 



MOND: iMilgroml Il983t iBegeman. Broeils. fc Sanders! 



[J99lt IMcGaugh fc de Blold 11998ft and so satisfy the 
first and second criteria. The third and fourth crite- 
ria are also satisfied in that the baryons in at least 
some outer disks have bee n shown to lie in a dynam- 
ical cold, undisturbed d isk (|Christlein fc Zaritskvi 120081 : 
iHerbert-Fort et aT1l2010ft . and the dark matter, while not 
in a perfectly spherical dist ribution, is thought to lie 
a significantly rounder one dHoekstra. Yee. fc Gladdersl 



l200l iMandelbaum et al.ll2006t ). 

For motions in the disk plane, the modified gravity 
and dark matter hypotheses are degenerate, with both 
predicting larger rotational velocities than inferred from 
the luminous matter using Newtonian gravity. The kine- 
matics in the vertical direction, where the baryons are 
strongly concentrated toward the disk plane, but the 
posited dark matter is not, provides the discrimination 
between t he two hypotheses. Such a test was imple- 
mented bv lSanchez-Salcedo et al.l (|2008l ) using the thick- 
ness of the H I disk at large radii in our own Galaxy. 
While they find a modest discrepancy with the MON- 
Dian prediction, they note that the comparison is some- 
what compromised by other potential sources of support 
for the gaseous layer, such as magnetic fields, cosmic 
rays, and SNe injected energy. Because one needs to find 
thicker than predicted H I to invalidate MOND, and be- 
cause uncertain sources of support will act to thicken the 
gas distribution, this test has an unresolved weakness. 
We propose to work in the same environment, which 
has so many appealing features, but to utilize another 
tracer population that cannot be supported by magnetic 
fields, cosmic rays, or injected energy. In this Letter 
we explore the implications of the measured dynamics 
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of outer disk star- forming knots flFerguson. et alj[l998t 
iThilker et aL|[200l iZaritskv fc Christleinll2007ft on mod- 
ified gravity models. 

2. TESTING THE MOND PARADIGM 

The vertical structure of galactic disks in modified 
Newtonian dynamics has been explored by a number of 
authors (e.g., Famaey & Binney 2005; Nipoti et al. 2007; 
Sanchez-Salcedo, Saha, & Narayan 2008). Using the phe- 
nomenological equation (Bekenstein & Milgrom 1984) 
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for the gravitational field $, where \i ( '—^ J is the phe- 

nomenological factor that depends on the magnitude of 
the local acceleration, g = |V$|, in units of the MOND 
scale ao — 10 -8 cm s -2 , Sanchez-Salcedo et al. (2008) 
showed that the disk scale-height in the deep MOND 
regime is given implicity by 
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In the above expressions, G is the gravitational constant, 
a is the velocity dispersion of matter in the vertical (out- 
of-disk plane) direction, El is the local column density of 
luminous (i.e., baryonic) material, v c is the circular ve- 
locity of matter, R is the galactocentric radius, and L(x) 
is the logarithmic derivative of fj, at x, where x = g/a^. 
Various fitting functions for /i(x), with the appropriate 
limiting behavior, have been developed, including what 
is often called the standard one, u — xj\J\ + x 2 (cf. 
iBeeeman. Broeils. fc Sanderslll99lD . and \x — x/(l + x) 
(|Famaev & Binnevil2005f )~~ 

We estimate the magnitude of the phenomenological 
factor \x at the radius of the knots using the equation for 
the radial acceleration 



/i (x) v 2 = 



GM L 
R 



(3) 



where v c is the circular velocity of matter, and Ml is 
the enclosed mass of luminous matter. This equation is 
valid for a spherical distribution of mass Ml. In fact, 
we are discussing a thin disk of baryonic material with a 
surface brightness (and corresponding mass column den- 
sity) that declines exponentially with radius, as is typical 
for spiral galaxies. For such a mass distribution the pre- 
dicted circular velocity is somewhat larger than that of 
a spherical distribution o f equivalent enclosed mass (cf. 
iBinnev k. Tremainelll987T ). which implies an even lower 
value of n and hence a stronger MOND effect. However, 
for large R the difference is small, < 10% at R > 5 disk 
radial scalengths. We confirm, using the standard fitting 
formula for n(x) that equation (3) produces a somewhat 
larger value of fi than the fitting formula (0.38 vs. 0.32, 
respectively) and therefore a correspondingly somewhat 
larger value of h — which is conservative given the argu- 
ment we are constructing. Combining equation (J2|) and 
$5§, we obtain 
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which we solve for h. 



The test is then a simple comparison of the predicted 
and observed values of h, which requires knowledge of the 
six measurable quantities: h, Ml, a, _R, E^, and v 2 . The 
more difficult of these to measure are h and a, and in par- 
ticular it is difficult to measure both quantities for any 
one galaxy because h is best measured in inclined sys- 
tems while a is best measured in face-on systems. One 
can combine characteristic values of h and/or <j to stand 
in place of the value for an individual galaxy in order 
to pursue this test. Of the two quantities, the one that 
has been measured in a number of systems is the s cale- 
height, h (jO'Brien. Freeman, fc van der Kruit]|2010[ ). al- 
though only for the H I layer and not yet for the knots 
themselves. This shortcoming must eventually be reme- 
died because otherwise our test is also susceptible to the 
unknown heating sources mentioned previously. Because 
the test has a stronger sensitivity to a, we present a mea- 
sure of a in the nearly face-on galaxy NGC 628 and use 
a representative estimate of h for disk galaxies. 

3. THE DATA 

We now describe how we estimate each of the quanti- 
ties necessary to evaluate equation ((4]) . 

3.1. The Disk Scaleheight, h 

The only currently available measurement of the scale- 
height of outer disks comes from measurem ents of the 
H I flaring (ISancisi k Al!eniri 979; M crrifieldlH992T: IQllind 
Il996t IQ'Brien. Freeman, fc van der Kruitl 1201(1 7 These 
of course (in the current context) have the sam e weak- 
ness described by iSanchez-Salcedo et al.1 (|2008l ) in that 
the thickness may be due to effects other than ther- 
m al pressure. Nevertheless, the currently defi nitive study 
of IQ'Brien. Freeman, fc van der Kruitl (|201Q[ ) shows that 
the typical scaleheight at the radii of interest here is ~ 
1 kpc. Their results also show a universal, steep rise in 
h toward the edge of the H I disk, suggesting that this 
test can become increasingly powerful at larger radii, or 
that additional heating is occurring near the edge of the 
neutral gas disk. Of the galaxies studied, IC 2531 is the 
nearest analog to NGC 628. They are both classified as 
Sc in NASA/IPAC Extragalactic Database. The V band 
magnitudes are —20.9 and —20.2, respectively. The H 
I extents are in both cases ~ 30 kpc and the rotation 
speeds are 220 and 170 km sec -1 . At the corresponding- 
mean radius of our observations (18 kpc; see below), h — 
1.2 kpc. The value of ft. as a function of radius for the two 
galaxies in their sample with circular velocity v c > 100 
km s -1 is the same, except for the final steep rise just 
before the end of the H I data. The uncertainty, judging 
from the scatter in h as a function of R and a comparison 
of the two more massive galaxies in the sample, appears 
quite small (~ 10%). The more serious, but not easily 
quantified, uncertainties relate to how well these galaxies 
describe NGC 628 and how well does the H I scaleheight 
describe that of the knots. With a typical velocity of 11 
km sec -1 , Ha-emitting knots will travel about 100 pc in 
their lifetime, significantly less than one scale height. 

The future of this test lies in directly measuring the 
scaleheight of the star forming knots. This can be done 
with either deep GALEX or Ha observations of edge- 
on galaxies, but has not yet been attempted. The value 
of the GALEX knots is that in their lifetime (> several 
hundred Myr) they will have completed several vertical 
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oscillations and will therefore test whether the H I and 
knots settle into the same vertical distribution. 

3.2. The Enclosed Mass of Luminous Matter, Ml 

The luminosity of the galaxy (M v — —20.73, based on 
m v = 9.25 (NED) and our adopted distance of 9.9 Mpc) 
is 1.6 x lO 10 ^© for M ViQ = 4.79. For B - V = 0.51 
and from the relations hips f or stellar mass-to-light ratios 
from lBell fc de Jond (|2001h . we estimate M/L v = 1.9 
and calculate that the stellar mass is 3.1 x 10 10 Mq. The 
uncertainty in this value will be dominated by that in 
M/L. The integrated magnitude is not likely to be off 
by more than 0.1 mag, and the scatter in distance mod- 
els (NED) corresponds to 0.1 mag as well. As a result 
the uncertainty in the luminosity cannot exceed 20%. 
Howev er, M/L could in princ iple be off by substantially 
more. IBell fc de Jond (|2001f) estimate, both from the 
differences obtained with different models and the scat- 
ter in their Tully-Fisher relationship, that the random 
scatter in logM/L is 0.2 (or about 50% in M/L, which 
is the uncertainty we adopt for Ml ) ■ To this mass w e 
add the gas as measured bv lKamphuis fc Briggsl (|1992j ). 
They measured a total H I mass of 12 x 10 10 M Q , of which 
they attribute 3 x 10 10 M Q to the distorted disk that lies 
beyond our measured knots. When we then correct the 
remainder for the presence of He by multiplying by 1.4, 
we estimate a gas mass of 12.6 x 10 10 Mq. The total 
mass we use is then 43.6 x 10 10 Af Q , with an uncertainty 
of 50%. 

3.3. The Velocity Dispersion, a 

The observational data consist of the radial velocity 
measurements of star-forming regions in the outer disk 
of the nearly face-on galaxy NGC 628. Star forming re- 
gions generate strong emission lines (Ha was used in this 
case) whose redshifts can be measured precisely. We cur - 
rently have a measurement by iHerbert-Fort et al.l (|2010l ) 
of the vertical velocity dispersion in the outer disk of 
NGC 628 using 14 such regions at galactocentric radii 
ranging from 1.04 to 1.79 R 25 (13 to 23 kpc). There 
are three principal sources of uncertainty in this mea- 
surement. First, the sample size is small. This could, 
in principle, be corrected and does not pose an inherent 
limit on the usefulness of this test. We have an esti- 
mate of thi s random uncertain t y , ^ 4 km sec -1 , from 



the work of IHerbert-Fort et al.l ( 201C0 . Second, there is 
an uncertainty regarding the orientation and kinematic 
parameters of the d isk about which we a re cal culating the 
dispersion. Again, IHerbert-Fort et all (|2010ft attempted 
to explore this uncertainty and concluded that for sev- 
eral of the more plausible orientations the uncertainty 
in a was about 2 km s . For the current purpose we 
adopt an uncertainty of 4 km s" 1 , but caution that there 
could be a systematic problem. This possibility should 
be addressed both with more data on this galaxy and 
with observations of other galaxies. 

A more fundamental limitation of the current data is 
that we measure the velocity dispersion only for knots 
that are relatively young (< 10 Myr) because we require 
the presence of the Ha emission line. This results in our 
using knots that are closely related to the neutral hydro- 
gen, and which have had insufficient time to erase the 
possible influence of unknown heating sources on the gas. 



The eventual solution to this problem is to measure the 
scaleheight of both this young population of knots and of 
the older GALEX population. If these two populations 
share a similar scaleheight, then the concerns of using 
a young tracer population are alleviated. Otherwise, we 
will need to somehow measure the velocity dispersion of 
the older, fainter knots. 

3.4. The Radius, R 

Due to the sparseness of the velocity field sampling, 
we must use data spanning a range of radii. The mean 
radius of the Herb ert-Fort et aLl (|2010( ) data is 18 kpc. 
The uncertainty in this comes from the uncertainty in 
the distance, which we estimate to be 10%. The more 
significant, and unquantifiable uncertainty comes from 
applying this test to data over a range of radii. In the 
future, a denser sampling of the kinematics will allow 
us to mitigate against potential systematic errors arising 
from this averaging (as well as enabling a stronger test 
by examining the radial behavior of any discrepancies in 
equation [4]). 

3.5. The Column Density of Luminous Matter, T,l 

Over this same region there are independent measure- 
ments of both the gaseous and stellar mass surface den- 
sities for NGC 62 cQ The H I gas, measure d via the 21 
cm line emission (|Kamphuis fc Br iggs 1992), has a mass 
of ~ 6 x 10 9 M Q over the relevant radial range, result- 
ing in an estimated mean surface gas mass density of 
5.6 Mq pc -2 after we multiply by 1.4 to account for 
He. This estimate of the gas mass excludes molecular 
gas, which we know is present at some level since there 
is ongoing star formation. The stellar surface density 
can be estimated from the measured surface brightness 
at 1.3i?25 (mv ~ 27 mag arcsec -2 ; iNatail et al.lll99lT ) 
and an estimated M/Ly of 3 for an old stellar popula- 
tion (for B — V = 1 and a modified Salpeter IMF, see 
IBell fc de Jo"ngl 12001 to be 1.0 M Q pc~ 2 . The surface 
mass density, £l, is therefore likely to be > 6.6 M Q pc -2 
at the representative radius of 1.3i?25- 

The uncertainty in this value is difficult to quantify in 
part because the data are not ours and in part because it 
is likely to be dominated by systematic errors. Because 
the H I mass dominates the sum, the principal source 
of uncertainty is likely to come from that measurement. 
We use the mean surface density evaluated over a large 
area. It would clearly be preferable to have more local 
measurements of both a and £l to avoid such averaging. 
We will assign the rather large uncertainty of 30% in an 
attempt to capture these possible errors. Ultimately, this 
test will require a more refined measurement of El and 
its uncertainty. 

3.6. The Rotational Velocity, v c 

The rotation curve is poorly measured because of the 
galaxy's low inclination. The rotational velocity at the 
outermost measured radius, ~ 12 kpc, is ~ 180 km s -1 
with no evidence for a decline with radius (jFathi et al.1 

1 The numerica l valu es presented here differ from those in 
IHerber t-Fort ct al. (2010) because we have 'conservatively' cho- 
sen to adopt a longer distance, which mitigates the discrepancy 
between the MONDian prediction and observations. 
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2007), although the velocity derived from their fitted 
mo del at the same radius is 15 km s _1 . The H I study 
by iKamphuis fc Briggs (|1992f ) adopted a flat rotation 
curve with amplitude 200 km s _1 and found good agree- 
ment with their data (except for a couple of high velocity 
clouds). We adopt v c = 170 km s" 1 with an uncertainty 
of 20%. 

4. COMPARISON TO PREDICTIONS 

Solving equation ^ using the discussed measured val- 
ues for NGC 628, and accounting for uncertainties using 
Monte-Carlo techniques, we find that h is predicted to be 
356^235 an d 298^172 kpc for the standard and Famaey- 
Binney interpolation formulae, respectively. The ob- 
served value of 1.2 kpc is rejected with slightly > 97% 
confidence in both cases. In a scenario with Newtonian 
gravity and in which the dark matter is not concentrated 
in the disk, so that the dominant mass component is 
the observed baryons, the predicted scaleheight, given 
by h = o- 2 /GvrI] L , is 1280^426 kpc, which is in agree- 
ment with the observations. 

There arc two ways in which we can improve the 
fidelity of this test. First, we can concentrate on 
lowering the uncertainties. Several of the terms con- 
tribute comparably to the final error, although the er- 
ror on a is currently the largest. It would be possi- 
ble to reduce this uncertainty by observing more Ha 
knots and using somewhat higher resolution observa- 
tions to minimize the importance of the instrumental 
contribution to the observed deviations. For compar- 
ison, t he H I velocity dispersions for the outer disks 
of the lO'Brien. Freeman, fc van der Kruitl (J2010I ) sam- 
ple are all ~ 7 km s _1 . These results suggest, but do 
not demonstrate, that the knot er's might be somewhat 
lower than our measurement. The next largest contribu- 
tion comes from the uncertainty in the rotation velocity. 
An analysis of the H I kinematics in which the rotation 
is solved for would help address this issue. Second, we 
can focus at larger radii, where the radial behavior differ- 
ences of the modified gravity and dark matter scenarios 



would provide additional discriminatory power. The dif- 
ficulty will be measuring a well in a small annulus since 
the Ha knots are somewhat rare. A realistic expectation 
is that we could reduce all of the uncertainties, except 
in the distance, by a factor of two given a feasible in- 
vestment of current-day observational resources. In that 
case, the discrepancy (for the current parameter values) 
would exceed the 99.99% confidence level. 

5. CONCLUSIONS 

We present the first use of the kinematics of outer disk 
star forming knots as a test of modifications to grav- 
ity at low accelerations. We evaluate the prediction 
for the disk scaleheight in m odified Newtonian gravity 
(jSanchez-Salcedo et al.l 120081 ) and compare to the typ- 
ical H I scaleheight of the outer disk in similar galax- 
ies. For our estimated uncertainties, we find a marginal 
(> 97% confidence) discrepancy between modified grav- 
ity scenarios and the observed scaleheight and excellent 
agreement between a model with Newtonian gravity and 
a baryonic disk. There remain various aspects of this test 
that require additional observations to confirm. In par- 
ticular measurements of the vertical scaleheight of the 
Ha-emitting knots and of the older GALEX-identificd 
knots, would directly address several potential system- 
atic uncertainties in the current version of the test. If 
those questions are satisfactorily resolved, we conclude 
that with additional data, this test has the potential to 
provide strong constraints (> 4er) on existing alternative 
gravity models with a relatively modest investment of 
observational resources. 
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